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Abstract 



A search for pair-produced charged Higgs bosons was performed in the high 
energy data collected by the DELPHI detector at LEP 11 at centre-of-mass 
energies from 189 GeV to 202 GeV. The three different final states, tvtv, 
cscs and cstv were considered. New methods were applied to reject wrong 
hadronic jet pairings and for the tau identification, where a discriminator based 
on tau polarisation and polar angles was used. No excess of data compared to 
the expected Standard Model processes was observed and the existence of a 
charged Higgs boson with mass lower than 71.5 GeV/c^ is excluded at the 95% 
confidence level. 
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1 Introduction 

The existence of a charged Higgs boson doublet is predicted by several extensions 
of the Standard Model. Pair-production of charged Higgs bosons occurs mainly via s- 
channel exchange of a photon or a Z° boson. In two-doublet models, the couplings are 
completely specified in terms of the electric charge and the weak mixing angle, 6\y, and 
therefore the production cross-section depends only on the charged Higgs boson mass. 
Higgs bosons couple to mass and therefore decay preferentially to heavy particles, but the 
details are model dependent. We assume that at LEP energies tUj. pair and a cs quark 
pair channels saturate the charged Higgs boson decays, and analyses of the three possible 
final states, tutu, cscs and csrz/, have been performed and are described in this paper. 
The Higgs decay branching fraction to leptons has been treated as a free parameter in 
the combination of the results of these three analyses. 

A search for pair-produced charged Higgs bosons was performed in the data collected 
by DELPHI during the LEP runs at centre-of-mass energies from 189 GeV to 202 GeV. 
The results reported here update those obtained in an earlier analysis of the DELPHI 
data limited to the 183 GeV run |jl|. Similar searches have been performed by the other 
LEP experiments [Q. 

A new technique was developed to improve the discrimination against the hadronic W 
decays in the search for H^ candidates. Improved methods using the r polarisation and 
boson production angles in the leptonic and semileptonic channels were used for rejection 
of W"'"W~ background. 



2 Data Analysis 

Data collected during the years 1998 and 1999 at centre-of-mass energies from 189 
GeV to 202 GeV were used. The total integrated luminosity of these data samples is 
approximately 380 pb~^. The DELPHI detector and its performance have already been 
described in detail elsewhere 

Signal samples were simulated using the HZHA generator P|. The background esti- 
mates from the different Standard Model processes were based on the following event 
generators: PYTHIA [| for gg(7), KORALZ [0] for and r+r", BABAMC [| for 

e^e~ and EXCALIBUR for four-fermion final states. Two-photon interactions were 
generated with TWOGAM |1^ for hadronic final states, BDK |11] for electron final states 



and BDKRC |]TT| for other leptonic final states. 

In all three analyses the final background rejection was performed by using a likelihood 
technique. For each of the N discriminating variables, the fractions F/^^^Xi) and F^^^{xi) 
of respectively H+H" and background events, corresponding to a given value Xi of the 
jth variable, were extracted from samples of simulated H'^'H" and background events 
normalised to equal size. The signal likelihood was computed as the normalised product 
of these individual fractions, [1^=1 n F!''' {xd I n Fl^'^ix,) + n.=i n F-''{x,)). 



2.1 The leptonic channel 

The signature for H+H^ — > t^Vt-t^Vt is large missing energy and momentum and two 
acoUinear and acoplanar jets containing either a lepton or one or a few hadrons. Tight 

'^The co-ordinate system used has the z-axis paraUel to the electron beam, and the polar angle calculated with respect 
to this axis. 

^The acoplanarity is defined as the complement of the angle between the two jets projected onto the plane perpendicular 
to the beam. 
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requirements for efficient operation of the most important sub-detectors were used in this 
analysis in order to ensure good quahty of the tracks. These result in slightly smaller 
integrated luminosities than in the hadronic channel (see Table ^). 



2.1.1 Event preselection 

To select leptonic events a total charged particle multiplicity between 2 and 6 was re- 
quired. All particles in the event were clustered into jets using the LUCLUS algorithm 
{djoin = 6.5 GeV/c) and only events with two reconstructed jets were retained. Both jets 
had to contain at least one charged particle and at least one jet had to contain not more 
than one charged particle. The angle between the two jets was required to be larger than 
30°. 

Two-fermion and two-photon events were rejected by requiring an acoplanarity larger 
than 13° if both jets were in the barrel region (43° < 9 < 137°) and larger than 25° 
otherwise. 

The two-photon background was further reduced by the following energy and momen- 
tum requirements: the sum of the jet energies multiplied by the sines of their jet angles 
to the beam direction, E^, was required to be larger than O.OSi/i if both jets were in the 
barrel region and larger than 0.1^^ in other cases; the total transverse momentum, 
to be greater than 0.04^^5; the total energy detected within 30° around the beam axis to 
be less than O.l^i; and the total energy outside this region to be greater than O.l^i. 

Additional r identification cuts were applied to reject WW events where the W's have 
not decayed to tv. If the r jet was identified as an electron it had to have a momentum 
below 0.13A/i and an electromagnetic energy below 0.14^^. For muons the momentum 
had to be below 0.13y^. If a r decay candidate particle was not identified as either a 
muon or an electron, it was considered to be a hadron and accepted as a r decay particle 
without further requirements. Events in which the invariant mass of either of the jets 
was more than 3 GeV/c^ were rejected. 

The effects of the r+z/^r^z/^ selection cuts are shown in Table for the combined 
189-202 GeV sample. 



2.1.2 Final background discrimination 

After these selections most of the remaining background consists of W"'"W~ 
T^v^T'Vr events. Events from both the H+H" signal and the W~*"W~ background have 
similar topologies and due to the presence of missing neutrinos in the decay of each of 
the bosons, it is not possible to reconstruct the boson mass. There are two important 
differences, however, that were used in order to discriminate the signal from the W"'"W~ 
background: the boson polar angle and the r polarisation. 

Assuming that the has a definite helicity, the polarisation (P^) of tau leptons orig- 
inating from heavy boson decays is determined entirely by the properties of weak in- 
teractions and the nature of the parent boson. The helicity configuration for the signal 
is H~ — i> Tj^i^TR (H"*" 'Tl^tl) for the W^ boson background it is W~ 'T'l^tr 
(W^ — > T^Vri) resulting in = +1 and = —1. The angular and momentum distri- 
butions depend on polarisation and it is possible to build estimators of the r polarisation 
to discriminate between the two contributions. 

The r decays were classified into the following categories: e, fi, n, n + nj, 3n and 
others. The information on the r polarisation was extracted from the observed kinematic 
distributions of the r decay products, i.e. their angles and momenta. These estimators 
are equivalent to those used at LEP I |T2[. For charged Higgs boson masses close to the 
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threshold, the boost of the bosons is relatively small and the r energies are similar to 
those of the r's from Z° decays (40-50 GeV). 

A likelihood to separate the signal from the W'^W^ background was built using four 
variables: the estimators of the r polarisation and the polar angle of the decay products 
of both r's. The distribution of that likelihood for data, expected backgrounds and a 
75 GeV/c^ charged Higgs boson is shown in Fig |I[ 



cut 


data 


total bkg. 


4-fermion 


other bkg. 


£^75 


Leptonic selection 


106274 


108275 


591 


107684 


73.6% 


Acoplanarity cut 


10841 


10519 


452 


10068 


60.0% 


Energy /momentum cuts 


360 


359 


343 


16 


45.6% 


r identification 


39 


39.7 


36.3 


3.4 


34.0% 



Table 1: The total number of events observed and expected backgrounds in the leptonic 
channel after the different cuts used in the analysis. The last column shows the efficiency 
for a charged Higgs boson signal with mH± = 75 GeV/c^. 



2.2 The hadronic channel 

In the fully hadronic decay channel, each charged Higgs boson is expected to decay 
into a cs pair, producing a four-jet final state. The two sources of background in this 
channel are the qqgg QCD background and fully hadronic four-fermion final states. In 
the four-fermion state background the significance of the Z^Z'^ pairs in the analysis is 
very small compared to the W^W~ pairs. This is due to the lower cross-section for Z^Z^ 
pair production and because the reconstructed Z°Z° pair masses are concentrated at high 
masses, around the Z° boson mass, out of the sensitivity reach of the analysis. Therefore, 
the four-fermion sample is referred to as W"''W~ in the rest of the paper. 



2.2.1 Event preselection 



Events were clustered into four jets using the Durham algorithm ||T3|. The particle 
quality requirements and the first level hadronic four-jet event selection followed in this 
analysis were the same as for the DELPHI neutral Higgs analysis ||14|| . 

In order to reject three-jet like QCD background events more effectively, the Durham 
clustering parameter value for transition from four to three jets (1/4^3) was required to 
be greater than 0.003. Events with a clear topology of more than four jets were rejected 
by requiring the 1/5^4 value for transition from five to four jets to be below 0.010 because 
of their worse di-jet mass resolution after forcing them into four jets. 

Energy-momentum conservation was imposed by performing a 4-C fit on these events 
and the difference between the two di-jet masses for each jet pairing was computed. A 
5-C fit, assuming equal boson masses, was applied in order to improve the di-jet mass 
resolution. The di-jet combination giving the smallest 5-C fit selected for the 

mass reconstruction. Events for which the 5-C fit divided by the number of degrees 
of freedom exceeded 1.5 or the difference of the masses computed with the same pairing 
after the 4-C fit exceeded 15 GeV/c^ were rejected. 
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2.2.2 Final background rejection 

The largest contribution to the part of the selected sample of W"'"W~ events whose 
reconstructed mass is below the W mass peak comes from picking one of the wrong di-jet 
pairings. These wrongly paired events are characterised by a larger difference between the 
masses of the two di-jets, i.e. the two boson candidates. As the initial quark antiquark 
pairs are connected by a QCD colour field, in which the hadrons are produced in the 
fragmentation process, the wrongly paired events can also be identified using a method 
of colour connection reconstruction |jl5| . 



The colour connection reconstruction method is based on the fact that, in the rest 
frame of the correctly paired initial quark antiquark pair, the hadrons that are produced 
in this colour string should have small transverse momenta relative to the quark antiquark 
pair axis. This could be distorted by hard gluon emission but such events are suppressed 
with the 1/5^4 Durham parameter cut. When boosted into a rest frame of a wrongly 
paired quark pair the transverse momenta of the particles relative to the quark quark 
axis are larger. The correct pairing is found by calculating the sum of transverse particle 
momenta in each of the three possible pairing hypotheses. The pairing chosen using the 
colour connection reconstruction is compared to the pairing chosen using the minimisation 
of the of the 5-C kinematical fit. The output of this comparison, called p_L-veto, is 
either agreement or disagreement and it is used later in the analysis as one of the variables 
in the background rejection likelihood. 

The production polar angle of the positively charged boson discriminates between 
W"''W~ and Higgs pairs. This angle is reconstructed as the polar angle of the di-jet with 
the higher sum of jet charges, where the jet charge is calculated as a momentum weighted 
sum of the charges of the particles in the jet The distribution of this variable allows 



the discrimination of the signal from the background of wrongly paired W"'"W~ events 
and QCD events, even though in these latter cases the variable does not correspond to a 
true boson production angle. 

Since the charged Higgs boson is expected to decay to cs in its hadronic decay mode, 
the QCD and W"^W~ backgrounds can be partially suppressed by selecting final states 
consistent with being cscs. A flavour tagging algorithm has been developed for the study 
of multiparton final states 0. This tagging is based on nine discriminating variables: three 
of them are related to the identified lepton and hadron content of the jet, two depend 
on kinematical variables and four on the reconstructed secondary decay structure. The 
finite lifetime of c (charm) particles is exploited to distinguish between c and light quark 
jets, while the c mass and decay multiplicity are used to discriminate against b jets. 
Furthermore s and c jets can be distinguished from u and d jets by the presence of 
an identified energetic kaon. Charged hadrons have been identified using the combined 
response of RICH and TPC dE/dx \TE\. The responses of the flavour tagging algorithm 



for the individual jets are further combined into an event cscs probability. 

The four variables described above: di-jet pair mass difference, the p^-veto, di-jet 
momentum polar angle and event cscs probability, were combined to form an event anti- 
WW likelihood function separating W"'"W~ events from H"'"H~ events. The response of 
this likelihood also discriminates H+H~ from the QCD background events. 

An anti-QCD likelihood was formed using all four variables which were used in the 
anti-WW likelihood and, in addition, two variables which can be used to separate QCD 
background from pair-produced bosons: the clustering algorithm parameter 2/4^3 and the 
event acoplanarity. 

'*A similar jet flavour tagging technique has been used in a determination of \Vcs\ at LEP II |^^. 
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The effects of the different sets of cuts are shown in Table ^ for the combined 189- 
202 GeV sample. The distribution of the anti-QCD likelihood on the preselection level 
and the distribution of the anti-WW likelihood after a cut on the anti-QCD likelihood 
are shown in Fig 0. The reconstructed mass distribution for data, expected backgrounds 
and signal after the anti-QCD and anti-WW cuts is shown in Fig ^ 



cut 


data 


total bkg. 


4-fermion 


other bkg. 


£^75 


4-jet presel. 


4156 


3981.8 


2515.0 


1466.8 


84.0% 


Durham y^ut 


3066 


2945.0 


2110.0 


835.0 


71.4% 




2397 


2340.0 


1760.5 


579.5 


60.5% 


Mass diff. 


1755 


1708.6 


1365.0 


343.6 


50.0% 


anti-QCD 


857 


844.5 


767.0 


77.5 


36.5% 


anti-WW 


653 


645.7 


577.8 


67.9 


33.2% 



Table 2: The total number of events observed and expected backgrounds in the hadronic 
channel after the different cuts used in the analysis. The last column shows the efficiency 
for a charged Higgs boson signal with mH± = 75 GeV/c^. 



2.3 The semileptonic channel 

In this channel one of the charged Higgs bosons decays into a cs quark pair, while the 
other decays into rVr- Such an event is characterised by two hadronic jets, a r candidate 
and missing energy carried by the neutrinos. The dominating background processes are 
QCD qqg event production and semileptonic decays of W"'"W~. The same requirements 
for efficient operation of the most important sub-detectors were used as in the analysis 
of the leptonic channel. 

2.3.1 Event preselection and r selection 

At least 15 particles, of which at least 8 were charged, were required. The total energy 
of the observed particles had to exceed 0.30A/i. The missing transverse momentum had 
to be greater than 0.08A/i and the modulus of the cosine of the angle between the missing 
momentum and the beam had to be less than 0.8. Events were also required to have no 
neutral particles with energy above 40 GeV. 

After clustering into three jets using the Durham algorithm, the clustering parameter 
?/3^2 was required to be greater than 0.003, and each jet had to contain at least one 
charged particle. The jet with the smallest charged particle multiplicity was treated as 
the r candidate and if two or more of the jets had the same number of charged particles, 
the jet with smallest energy was chosen. The r candidate was required to have no more 
than six particles, of which no more than three were charged. 

2.3.2 Final background rejection 

The mass of the decaying bosons was reconstructed using a constrained fit requiring 
energy and momentum conservation with the known beam energy and imposing the 
masses of the two bosons to be equal. The three components of the momentum vector of 
the Vr and the magnitude of the r momentum were treated as free parameters, reducing 
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the number of degrees of freedom of the fit from 5 to 1. Only events with a reconstructed 
mass above 40 GeV/c^ and a below 2 were selected. 

Separate likelihood functions were defined to distinguish the signal events from the 
QCD and the W"'"W~ backgrounds, in a manner similar to that used for the other channels 
described above. 

To define the anti-QCD likelihood, the acollinearity of the event (after forcing the 
event into two jets), the polar angle of the missing momentum, the logarithm of the 
clustering parameter 1/3^2, and the product of the r jet energy and the smaller of the two 
angles between the r jet and one of the other jets were used as discriminating variables. 

For the event anti-WW likelihood the variables used were the reconstructed polar 
angle of the negatively charged boson (where the charge was determined from the leading 
charged particle of the r jet), the angle between the boson and the r in the boson rest 
frame, the energy of the r jet, the classification of the decay of the r candidate (e, /i, vr, 
TT + n7, Stt and others), and the cs probability of the hadronic di-jet. 

The effects of the different sets of cuts are shown in Table ^ for the combined 189- 
202 GeV sample. The distribution of the anti-QCD likelihood on the r selection level 
and the distribution of the anti-WW likelihood after a cut on the anti-QCD likelihood 
are shown in Fig |. The reconstructed mass distribution for data, expected backgrounds 
and signal after anti-QCD and anti-WW cuts is shown in Fig ^. 



cut 


data 


total bkg. 


4-fermion 


other bkg. 


£^75 


preselection 


6395 


6104.8 


3043.8 


3061.0 


81.6% 


r selection 


2149 


2144.0 


1788.0 


356.0 


58.1% 




1667 


1699.0 


1552.4 


146.6 


50.4% 


likelihoods 


325 


307.4 


292.4 


15.0 


33.8% 



Table 3: The total number of events observed and expected backgrounds in the 
semi-leptonic channel after the different cuts used in the analysis. The last column 
shows the efficiency for a charged Higgs boson signal with mH± = 75 GeV/c^. 



3 Results 

3.1 Selection efficiencies and uncertainties 

The number of real data and background events and the estimated efficiencies for these 
selections for two different masses are summarised in Table ^ for the three final states. 
The quoted errors include the systematic uncertainties in the expected background and 
the signal efficiency. Small contributions to these uncertainties are due to uncertainties 
in the luminosity measurement and in the cross-section estimates of the generated Monte 
Carlo samples. 

The event selection and systematic errors in the leptonic analysis are very similar to 
those in the DELPHI leptonic W+W" analysis [|r^. The largest part of the background 
and signal efficiency uncertainties in the leptonic channel is due to the limited simulation 
statistics available. Combining these uncertainties gives a total uncertainty of the order 
of 10% in the background rate and 5% in the signal efficiency. 

The largest contribution in the hadronic and semileptonic analyses is due to differences 
in the distributions of the preselection and likelihood variables in data and simulation. 
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Chan. 




lum. 


data 


total bkg. 


^70 


S75 


TVTV 


189 


153.8 


16 


15.0± 1.5 


32.3±1.6% 


34.2±1.6% 


TVTV 


192 


24.5 


3 


2.8± 0.3 


33.6±1.6% 


34.2±1.6% 


TVTV 


196 


72.4 


10 


8.6± 0.8 


33.6±1.6% 


34.2±1.6% 


TVTV 


200 


81.8 


8 


9.0± 0.9 


33.6±1.6% 


33.5±1.6% 


TVTV 


202 


39.4 


2 


4.4± 0.4 


33.6±1.6% 


33.5±1.6% 


CSCS 


189 


154.3 


288 


267.8±16.1 


36.2±2.0% 


33.1±2.0% 


CSCS 


192 


25.5 


36 


42.7± 2.6 


36.2±2.0% 


33.1±2.0% 


CSCS 


196 


77.1 


141 


130.0± 7.8 


34.5±2.0% 


33.8±2.0% 


CSCS 


200 


83.9 


133 


138.8± 8.3 


33.9±2.0% 


33.5±2.0% 


CSCS 


202 


40.6 


55 


66.2± 4.1 


33.9±2.0% 


33.5±2.0% 


CSTV 


189 


153.8 


126 


118.8± 6.9 


33. Oil. 8% 


31.8±1.7% 


CSTV 


192 


24.5 


29 


21. 2± 1.9 


36.6±2.2% 


35.2±2.1% 


CSTV 


196 


72.4 


76 


62.2± 5.5 


36.6±2.2% 


35.2±2.1% 


CSTV 


200 


81.8 


67 


71. 1± 6.4 


35.2±2.1% 


35.3±2.1% 


CSTV 


202 


39.4 


27 


34.1± 3.1 


35.2±2.1% 


35.3±2.1% 



Table 4: Integrated luminosity, observed number of events, expected number of back- 
ground events and signal efficiency (70 GeV/c^ and 75 GeV/c^ masses) for different 
decay channels and centre-of-mass energies. 

The systematic error on the efficiency of the common DELPHI hadronic four-jet preselec- 
tion has been estimated to be ±4% HT^l. The uncertainties related to the other selection 



variables have been estimated by comparing the shapes of the variable distributions in 
data and simulation. This has been done at the preselection level where the background 
event rate is so large that a possible signal would have no effect on the global shapes of 
the variables. The agreement of all variables has been found to be satisfactory to the level 
of a few percent. Combining these errors, a total uncertainty of 6% has been estimated 
for the background rate and signal efficiency in the hadronic channel. In the semileptonic 
channel the combined background error estimate is 6 to 9% depending on the energy 
sample and the error of the signal efficiency is of the order of 6%. The combined error 
estimates are included in Table |[ 

3.2 Determination of the mass limit 

No significant signal-like excess of events was observed in any of the three final states 
investigated. We find an agreement between data and background expectations also in 
the mass region around 68 GeV/c^ where the L3 collaboration has reported an excess 
using data collected at the same centre-of-mass energies [Q]. A lower limit for a charged 
Higgs boson mass was derived at 95% confidence level as a function of the leptonic Higgs 
decay branching ratio BR(H — > tv^). The confidence in the signal hypothesis, CLg, was 
calculated using a likehhood ratio technique . 



The background and signal probability density functions of one or two discriminat- 
ing variables in each channel were used. The data samples collected at the five different 
centre-of-mass energies were treated separately in the combination. Correlations between 
the systematic uncertainties in different centre-of-mass energies were not included. In 
the hadronic and semileptonic channels the two discriminating variables were the recon- 
structed mass and the anti-WW likelihood; in the leptonic channel only one background 
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discrimination likelihood was used since mass reconstruction is not possible. The dis- 
tributions of the discriminating variable for signal events, obtained by the simulation at 
different mass values for each ^/s, were interpolated for intermediate mass values. To 
obtain the expected signal rate at any given mass the signal efficiencies were fitted with 
polynomial functions. 

A Gaussian smearing of the central values of the number of expected background 
events by their estimated uncertainties was introduced in the limit derivation program. 

The results are summarised in Fig ^. A lower mass limit of Mh± > 71.5 GeV/c^ 
can be set at the 95% confidence level, independently of the branching ratio BR(H 
TUr). The production cross-section for the charged Higgs bosons signal at this mass 
(71.5 GeV/c^) and a centre-of-mass energy of 202 GeV would be 0.255 pb. The median 
of the limits obtained from a large number of simulated experiments is 73.3 GeV/c^. 



4 Conclusion 

A search for pair-produced charged Higgs bosons was performed using the full statis- 
tics collected by DELPHI at LEP at centre-of-mass energies from 189 GeV to 202 GeV 
analysing the tutu, cscs and cstv final states. No significant excess of candidates was 
observed and a lower limit on the charged Higgs mass of 71.5 GeV/c^ is set at 95% con- 
fidence level. The sensitivity of the DELPHI charged Higgs boson search is comparable 



to the ones of the other LEP collaborations 21 
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Figure 1: Distribution of the anti-WW likelihood for leptonic events at 189-202 GeV. 
The expected histogram for a 75 GeV/c^ charged Higgs boson signal has been normalised 
to the production cross-section and 100% leptonic branching ratio and added to the 
backgrounds. 
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Figure 2: Distributions of the anti-QCD and anti- WW likelihoods for hadronic events 
at 189-202 GeV. The anti-QCD likelihood is plotted on the preselection level and the 
anti- WW likelihood after a cut on the anti-QCD likelihood. The generated H"'"H~ signal 
mass is 75 GeV/c^ and the signal histograms have been normalised to the production 
cross-section and 100% hadronic branching ratio, multiplied by a factor of 15 and super- 
imposed on the background histograms. The arrows indicate the cut values, below which 
events were rejected. 
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Figure 3: Reconstructed mass distribution of hadronic events at 189-202 GeV at the final 
selection level. The generated H"'"H~ signal mass is 75 GeV/c^ and the signal histogram 
has been normalised to the production cross-section and 100% hadronic branching ratio 
and added to the backgrounds. 
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Figure 4: Distributions of the anti-QCD and anti- WW likelihoods for semileptonic events 
at 189-202 GeV. The anti-QCD likelihood is plotted on the r selection level and the 
anti- WW likelihood after a cut on the anti-QCD likelihood. The generated H"'"H~ signal 
mass is 75 GeV/c^ and the signal histograms have been normalised to the production 
cross-section and 50% leptonic branching ratio, multiplied by a factor of 50 and superim- 
posed on the background histograms. The arrows indicate the cut values, below which 
events were rejected. 
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Figure 5: Reconstructed mass distribution of semileptonic events at 189-202 GeV at the 
final selection level. The generated H^H^ signal mass is 75 GeV/c^ and the signal his- 
togram has been normalised to the production cross-section and 50% leptonic branching 
ratio and added to the backgrounds. 
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Figure 6: The 95% confidence level observed and expected exclusion regions for H''' in 
the plane BR(if —>■ TV-r) vs. Mh± obtained from a combination of the search results 
in the fully leptonic, hadronic and semileptonic decay channels at \/s = 189-202 GeV. 
The expected median of the lower mass limits has been obtained from a large number 
of simulated experiments. The median is the value which has 50% of the limits of the 
simulated experiments below it and the ± Icr lines correspond similarly to 84% and 16% 
of the simulated experiments. 



